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Abstract 

The pseudohyphal growth response is a dramatic morphological transition and presumed foraging mechanism wherein 
yeast cells form invasive and surface-spread multicellular filaments. Pseudohyphal growth has been studied extensively as a 
model of conserved signaling pathways controlling stress responses, cell morphogenesis, and fungal virulence in 
pathogenic fungi. The genetic contribution to pseudohyphal growth is extensive, with at least 500 genes required for 
filamentation; as such, pseudohyphal growth is a complex trait, and linkage analysis is a classical means to dissect the 
genetic basis of a complex phenotype. Here, we implemented linkage analysis by crossing each of two filamentous strains 
of Saccharomyces cerevisiae (21278b and SKI) with an S288C-derived non-filamentous strain. We then assayed meiotic 
progeny for filamentation and mapped allelic linkage in pooled segregants by whole-genome sequencing. This analysis 
identified linkage in a cohort of genes, including the negative regulator SFLl, which we find contains a premature stop 
codon in the invasive SKI background. The S288C allele of the polarity gene PEA2, encoding Leu409 rather than Met, is 
linked with non-invasion. In El 278b, the peo2-M409L mutation results in decreased invasive filamentation and elongation, 
diminished activity of a Ksslp MARK pathway reporter, decreased unipolar budding, and diminished binding of the 
polarisome protein Spa2p. Variation between SKI and S288C in the mitochondrial inner membrane protein Mdm32p at 
residues 182 and 262 impacts invasive growth and mitochondrial network structure. Collectively, this work identifies new 
determinants of pseudohyphal growth, while highlighting the coevolution of protein complexes and organelle structures 
within a given genome in specifying complex phenotypes. 
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Introduction 

The budding yeast Saccharomyces cerevisiae undergoes a 
pronounced growth transition in response to nitrogen limitation or 
glucose limitation, forming multicellular pseudohyphal filaments that 
can spread outward from a colony and/ or invade the surface of a 
solid growth substrate [1,2]. Yeast pseudohyphal filament formation 
is a presumed foraging mechanism, accomplished through under- 
lying changes in cell adhesion, cell cycle progression, and budding 
[1,3,4]. During pseudohyphal growth, yeast cells remain physically 
connected after cytokinesis via mechanisms encompassing the 
regulated expression and shedding of the flocculin Flo lip [5-7]. 
Cells undergoing pseudohyphal growth exhibit increased apical 
growth through reorganization of the actin cytoskeleton, regulation 
of polarity proteins, and delayed G2/M progression [8-12]. 

The molecular basis of yeast pseudohyphal growth has been studied 
extensively as a model of conserved signaling pathways controlling cell 
morphogenesis and polarity. Furthermore, related processes of 
filamentous development in the principal opporainistic human fungal 
pathogen Candida albicans are required for virulence, and signaling 
pathways between the related yeasts are conserved [13]. Classic 



studies of pseudohyphal growth in S. cerevisiae have resulted most 
prominently in the identification of core pseudohyphal growth 
signaling modules encompassing the Ksslp mitogen-activated protein 
kinase (MAPK) cascade, the cAMP-dependent protein kinase A 
(PKA) pathway, and the AMP-activated protein kinase ortholog Snflp 
[14—20]. The pseudohyphal growth MAPK cascade encompasses 
Stel Ip, Ste7p, and the MAPK Ksslp [10,14]. Ksslp phosphorylates 
the Ste 1 2p transcription factor, resulting in dissociation of the negative 
regulatory Diglp and Dig2p interactors and binding of a Stel2p- 
Teclp heterodimer to target promoters, such as die FLOll promoter 
[21-23]. Tpk2p, a catalytic subunit of PKA, phosphorylates the Flo8p 
transcription factor, promoting Flo8p binding and transcriptional 
activation at the FLOll promoter and other regulatory sites [17,24— 
26]. In response to glucose limitation, FLOll transcription is 
regulated by Snflp; the Snflp-Gal83p isoform promotes cell adhesion 
during invasive filamentation by antagonizing Nrglp- and Nrg2p- 
mediated repression oi FLOll [19,27]. 

While the central components of these signaling pathways have 
been identified, the scope of the yeast pseudohyphal stress 
response is broad [28-33], and the mechanisms enabling these 
genes and gene products to drive pseudohyphal filamentation are 
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Author Summary 

Cellular processes in eukaryotes are brought about through 
the contributions of large gene sets, and a continuing 
obstacle in studying these processes lies in the identification 
of critical constituent genes. The yeast pseudohyphal growth 
transition is an Important example of a complex cellular 
growth transition. During pseudohyphal growth, yeast cells 
form connected chains or filaments, constituting a means of 
foraging for nutrients under conditions of nitrogen and/or 
glucose limitation. Yeast pseudohyphal growth has been 
studied for over two decades as a model of signaling systems 
controlling stress responses, cell shape, and fungal virulence. 
Hundreds of genes are required for pseudohyphal growth, 
however, and the critical genes that determine the filamen- 
tous phenotype have not been elucidated. Towards this goal, 
we implemented a genetic approach to identify alleles linked 
with the pseudohyphal growth phenotype. These studies 
identified previously unstudied variation in proteins func- 
tioning in a complex that controls cell polarity and in a 
protein of the mitochondrial Inner membrane. This work 
indicates that proteins in complexes and organelles have 
coevolved within a given genome to yield distinct outputs 
and phenotype, while highlighting the application of an 
approach that is useful for the analysis of complex 
phenotypes In many eukaryotes. 

incompletely defined, as are the genetic determinants within this 
gene set that underlie filamentation. To further dissect pseudohy- 
phal growth pathways, we implemented a linkage study, coupling 
whole genome sequencing with pooled segregant analysis. The 
results present previously unidentified genetic determinants of 
yeast invasive growth and indicate the coevolution of proteins 
within complexes in driving phenotype. 

Results 

Pooled Segregant Analysis of Yeast Invasive 
Filamentation with Deep Sequencing 

For linkage analysis, we selected as parents the non-filamentous 
S288C-derived strain BY4741 and the filamentation-competent 
strains SI 278b and SKI [34,35]. Filamentous-form growth in 
haploid strains is classically assessed using the plate-washing assay of 
Gimeno et al. [1] to identify pseudohyphal cells that have invaded 
the agar substrate. The invasive phenotype of each parent strain in 
this assay is indicated in Fig. 1 A. The experimental design of the 
linkage study is presented in Fig. IB. The non-invasive S288C- 
derived strain was mated with each of the filamentous SI 278b and 
SKI strains, and the resulting diploid strain from each cross was 
sporulated. Meiotic progeny from dissected tetrads were assayed for 
agar invasion by plate-washing, and spores indicating strongly non- 
invasive or invasive phenotypes were pooled for subsequent linkage 
analysis. Only spores resulting from complete meiosis were included 
in these phenotypic pools, and intermediate filamentation pheno- 
types were excluded from subsequent analysis to provide the 
greatest likelihood of identifying allelic variation with a strong effect 
on filament formation. Genomic DNA was extracted from each 
segregant pool and subjected to high-throughput sequencing that 
yielded greater than 100-fold coverage per pool. 

Identification of Allelic Variants Linked with Invasive 
Phenotype between a Non-filamentous S288C Derivative 
and SI 278b 

From the BY4741-by-S1278b cross, 31 complete tetrads (124 
spores) were screened for agar invasion, identifying 37 strongly 



invasive spores and 63 non-invasive spores (Fig. 2A). The 
segregant pools were sequenced, and candidate determinants of 
the invasive phenotype were identified using a linkage LOD score 
of greater than 3 as an arbitrarily defined cut-oflF. Table SI 
provides a listing of these alleles, encompassing only variants that 
are in protein-coding sequence and that are non-synonymous with 
respect to the encoded amino acid sequence. This allele set affects 
50 genes in eleven linkage blocks physically located on seven yeast 
chromosomes. Figure SI summarizes the available functional 
information for this gene set. 

Representative plots of non-synonymous allelic variation with 
respective LOD scores are graphed in Fig. 2B for chromosomes V 
and IX, highlighting the pseudohyphal growth transcription factor 
gene FLOS and the flocculin effector gene FLOll. FLOS is a 
pseudogene in S288C-derived strains [36], and in this analysis, the 
BY4741 allelic variant containing a premature translational stop at 
codon 142 of the FLOS sequence yielded a LOD score greater 
than 17 (Fig. 2B and C). The FLOll locus exhibits fifteen allelic 
changes linked with invasive growth phenotypes (Fig. 2B and C). 
Previous studies identified allelic variation in FLOll sequence 
encoding amino- and carboxy-terminal regions linked with the 
ability to form biofUms on the surface of wine [37]. We recovered 
these as well as additional sites of DNA sequence variation in 
FLOll, with the S1278b-encoded alleles indicating linkage with 
strong invasive growth. The FLOll sequence contains an internal 
repeat region that is a source of allelic variation between some 
strains and colonies [38,39]; however, we did not observe a change 
in the number of these repeats between BY4741 and SI 278b. 
Collectively, the identification of these known pseudohyphal 
growth genes demonstrates the relevance of results obtained from 
our pooled segregant analysis. 

PEA2 Allelic Variation Impacts Cell Morphology and 
MAPK Pathway Activity 

To further identify important determinants of invasion, we 
screened candidates from Table SI as follows: 1) we generated 
gene deletions and assayed for invasive growth phenotypes (Table 
S2), and 2) for genes yielding deletion phenotypes, we generated 
mutants with swapped alleles to identify genetic variants required 
for invasive growth in SI 278b. In particular, we focused on alleles 
of genes that contributed to cell polarity, cell cycle progression, cell 
morphology, and cell responses to nitrogen/carbon limitation, as 
these are hallmark characteristics of filamentation. 

By this approach, we identified variation in PFA2 as an 
important part of the genetics underpinning invasive growth. 
Pea2p localizes to sites of polarized growth as a component of a 
protein complex, termed the polarisome [40,41]. PEA2 is required 
for wild-type invasive growth, mating projection formation, and 
bipolar bud site selection in diploids [42,43]. In the filamentous 
SI 278b strain, PEA2 codon 409 specifies methionine rather than 
the leucine residue encoded in the S288C-derived reference 
genome. The pea2-M409 allele was linked with invasive growth, 
and generation of an integrated site-specific mutation (pea2- 
M409L) reconstituting the S288C-encoded allele in SI 278b 
resulted in decreased invasive growth (Fig. 3A). Relative to wild 
type SI 278b, the cell morphology of the pea2-MW9L mutant is 
altered, exhibiting decreased elongation (Fig. 3B); over a popula- 
tion of 200 cells, the percentage of ^«a2-M409L cells with a 
length:width ratio of less than 1.5 was nearly four-fold the 
corresponding percentage in a wild type strain. In addition, the 
^m2-M409L mutant is impaired in Ksslp MAPK signaling 
activity. The Ksslp kinase activates the Stel2p/Teclp transcrip- 
tion factor complex, which recognizes a regulatory element (FRE) 
in the FLOll promoter. The plasmid-based Pfio^^.g/io-iocZ 
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Figure 1. Overview of pooled segregant analysis to identify linkage with invasive growth phenotypes. A) Invasive growth phenotypes 
of the non-filamentous S288C-derived BY4741 strain and the filamentous Z1278b and SKI strains. Spotted cultures of the strains are shown on YPD 
medium before and after washing the surface of the agar plate under a gentle stream of water. Residual cells after washing are indicative of invasive 
growth. The degree of invasive growth is shown overlaid on each respective image and in the boxed inset, with "-H-" indicating strong invasive 
growth, and "-" indicating an absence of invasive growth. B) Schematic depicting the pooled segregant whole-genome sequencing approach. 
doi:1 0.1 371 /journal.pgen.1 004570.g001 



construct contains the Stel2p/Teclp-responsive region of ttie 
FLO 11 promoter fused to lacZ [6], and, by this reporter, the pea2- 
M409L mutant yields significantly decreased Stel2p/Teclp- 
dependent transcriptional activation of FLO 11 relative to wild- 
type Z 1278b (Fig. 3C). In contrast, the pea2-M409L mutation 
results in wild-type levels of a similarly designed FLO 11 promoter 
fusion responsive to the PKA pathway effector Flo8p (Fig. 3C) [6] . 

Genetic Variation in PEA2 IVlodulates Bud Site Selection 
and Spa2p Binding 

Under conditions of vegetative growth haploid yeast cells bud in 
an axial pattern, with new buds emerging adjacent to the preceding 
bud site [44]. Haploid cells undergoing pseudohyphal growth, 
however, adopt a predominantly unipolar budding pattern wherein 
the first bud forms distal to the original cell division site, and 
subsequent buds cluster in the distal pole [1,10]. Here, we find that 
in the SI 278b background the ^m2-M409L mutant, corresponding 
to the S288C-encoded PEA2 allele, exhibited a decrease in unipolar 
budding and an increase in axial budding relative to wild type (p< 
0.001), with levels intermediate between an otherwise isogenic wild- 
type strain and a pea2A mutant (Fig. 3D). For this analysis, we 
examined a population of invasive cells exhibiting three or more 
bud scars, such that patterns of axial, unipolar, bipolar, and random 
budding could be reliably distinguished [44,45]. This budding 
phenotype was evident in invasive cells, but not in an equally sized 
population of cells scraped from the surface of an agar plate. 
Previous studies have indicated that the majority of bud sites are 
distal in a pea2A mutant [12]; results here also indicate that the 
majority of bud sites are distal in pea2 mutants, but budding pattern 
analysis does indicate that Pea2p residue 409 impacts unipolar 
budding in invasive haploid cells. 

hi the polarisome complex, Pea2p binds the scaffolding protein 
Spa2p, a large coiled-coU domain-containing protein required for 
polarisome function [41,46]. Here, we assessed the possibility that 
allelic variation at the PEA2 locus impacts Spa2p binding, using 
Protein A (ProA)-tagged Pea2p variants to recover by co-immuno- 
precipitation Spa2p tagged at its amino terminus with the hemagglu- 
tinin (HA) epitope. In the Z 1278b strain, the Pea2p-M409-ProA 
variant recovered more HA-Spa2p than the Pea2p-L409-ProA 
variant, representing the S288C-encoded PFA2 allele (Fig. 3E). 

Analysis of Genetic Variation Linked with Invasion 
between BY4741 and SKI 

The BY4741-by-SKl cross was implemented as described in 
Experimental Procedures, and phenotypic analysis of meiotic 



progeny identified 51 and 24 strongly invasive and non-invasive 
spores, respectively (Fig. 4A). Subsequent deep-sequencing of the 
phenotypic pools identified allelic variation linked with invasive 
growth phenotypes within eleven separated locus blocks encom- 
passing 88 genes exhibiting non-synonymous changes and a LOD 
score of greater than 4 (Table S3). A functional breakdown of 
these genes is indicated in Figure S2. In this analysis, we used a 
higher LOD score relative to the 21278b-by-BY4741 cross in 
order to limit the number of selected allelic variants to a 
manageable size for further study, as the SKI and BY4741 
genomes are more divergent (99.5% sequence identity) than the 
21278b and BY4741 genomes (99.7% identity) [47]. Very few 
allelic variants linked with invasive growth in SI 278b were also 
identified in SK 1 , aside from a few sequences near FLOS that are 
unlikely to be causative. The set of identified alleles was primarily 
distinct between the two linkage studies, and deletion phenotypes 
for tested genes in SKI are indicated in Table S4. 

In haploid spores from this cross, genetic variation is most 
strongly linked with the invasive growth phenotype over a region of 
roughly 80,000 bp on chromosome XV, encompassing SFLl 
(Fig. 4B and C). SFLl encodes a transcriptional repressor of 
pseudohyphal growth that functions by binding to the FLO 11 
promoter, thereby blocking transcriptional activation [24,25,48- 
51]. Consistent with its function in repressing i^LO gene expression, 
deletion of SFLl results in exaggerated invasive growth [28,32,48]. 
Interestingly, the SKI strain contains an allelic variant ofSFLl with 
respect to S288C-derived strains, resulting in the conversion of 
codon 477 (CAA encoding glutamine) to a TAA stop codon 
(Fig. 4D). This premature stop codon truncates SFLl prior to the 
sequence encoding a domain (AA 571-658) that is strongly similar 
to a conserved region in Myc oncoproteins [48]. Previous studies 
have identified hyperactive fdamentation in a mutant of the 
CEN.PK 1 1 3-7D background upon introduction of a premature 
translational stop at SFLl codon 320 (Q,320-stop) [51]. Here, we 
found that allelic variation in SFLl , encompassing a premature stop 
codon (C 1430T, Q477-stop) in the SKI background, is linked to the 
aggressively invasive phenotype of SKI relative to BY4741. 

Allelic Variation in MDM32 Affects Mitochondrial 
Morphology and Invasive Growth 

In the BY4741-by-SKl cross, allelic variation in MDM32 was 
linked to invasive growth more strongly than any other identified 
locus, with a LOD score of 9. As indicated in Fig. 4B-D, MDM32 
is found on chromosome XV, and relative to BY4741, the SKI 
allele oi MDM32 encodes Serl82 and Phe262 rather than Cys 
and Leu, respectively. MDM32 encodes a protein complex 
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Figure 2. Identification of alleles linked with invasive growth in meiotic progeny from the BY4741-by-D1278b cross. A) Summary of 
phenotypic screening for invasive growth in haploid spores. B) Plots of non-synonymous allelic changes graphed by chromosomal location are 
presented, with allelic variation exhibiting a LOD score of greater than 3 highlighted in red. Clusters of alleles with LOD scores greater than 3 are 
boxed, and expanded blow-ups of loci on Chromosome V and IX housing the invasive growth determinants FLOS, PEA2, and FLOll are shown. Red 
arrows represent genes, with forward directionality indicating the Watson strand of the respective chromosome. LOD scores are plotted above each 
gene. The width of the LOD score bar for each gene corresponds to the chromosomal boundaries of the allelic changes observed in the respective 
gene coding sequence; a single nucleotide change between the sequences is indicated by a thin line. Chromosomal nucleotides delineating the 
boundaries of the highlighted region are shown at the upper right and left of each expanded plot. C) Specific nucleotide and amino acid changes 
between alleles of the FLOS, PEA2, and FLOll genes are indicated, along with the respective LOD score for each gene. 
doi:1 0.1 371/journal.pgen.1 004570.g002 



subunit of the mitochondrial inner membrane required for 
membrane organization, the maintenance of elongated mitochon- 
drial morphology, and mitochondrial DNA nucleoid stabilization 
[52]. Mitochondrial function is required for pseudohyphal growth, 
as filamentation-competent strains of S. cerevisiae containing a 
deleted version of the mitochondrial genome are unable to form 
pseudohyphae [53,54]; however, a role for Mdm32p in enabling 
invasive growth remains to be identified. 

To determine the effect oi MDM32 allelic variation on pseudohy- 
phal growtli, we replaced the SKl-encoded MDM32-C546/T787 



allele, specifying Mdm32p-S182/F262, with BY4741-encoded 
MZ)M52-G546/A787, specifying Mdm32p-C182/L262, in the SKI 
genetic background. This allelic swap decreased invasive growth in 
SKI, and agar invasion was rescued upon reintroduction of the native 
SKl-encoded MDM32 allele (Fig. 5A). The SKI mutant containing 
tlie BY4741 -encoded allele oi MDM32 exhibited a more rounded cell 
morphology, witli the percentage of cells displaying a cell lengtli: width 
ratio of less than 1.5 increasing from 17% in wild-type SKI cells to 
68% in the SKI mutant (Fig. 5B). Reintroduction of SKl-encoded 
MDM32 recovered levels of cell elongation similar to wild type. 
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Figure 3. Variation at Pea2p residue 409 impacts invasive growth, budding, and polarisome interactions. A) Invasive growth of the 
wild-type El 278b strain, pea2A strain, and a strain containing the integrated pea2-M409L point mutation in the El 278b background was assayed by 
plate washing (++, strong invasive growth; +, invasive growth; — , lack of invasive growth). Invasive growth levels were quantified by determining the 
ratio of pixel intensity of images of the spotted cultures before and after plate washing; these ratios are indicated to the lower right of each spotted 
culture image post-washing. B) Cell morphology was assessed by determining the length-to-width ratio of 200 cells each for the indicated strains. 
Prior to imaging, cells grown on YPD medium were scraped into a suspension in liquid media. For length-width calculations, the longest cell axis was 
selected for cell length measurements. C) Activity of MARK and PKA signaling pathways in pea2 mutants. The bar graphs indicate activity of plasmid- 
based cAMP-dependent PKA and pseudohyphal growth MAPK pathway-responsive Pfio,,-/ocZ reporters. The PKA pathway-responsive Pfici reporter 
contains a fragment of the FL011 promoter (Pfion-e/?) 21.0-21.4 kb from the FL011 start codon, and the Kssip MAPK-responsive Pfion reporter 
contains a promoter fragment (Pfioi 1-9/10) 21.6-22.0 kb upstream of the FLOll start codon. )3-galactosidase activity was measured from cells growing 
exponentially in SC-Ura medium. Error bars indicate standard deviations of three independent measurements. D) Budding patterns were analyzed for 
the indicated strains by calcofluor white staining of invasive cells recovered from the agar. Sample images of the observed budding patterns are 
presented. Budding patterns of the pea2A and pea2-M409L mutants are distinct from wild type, with p-values of <0.01 and <0.001, respectively, by 
Pearson's chi-squared test. E) Western blots indicate co-immunoprecipitation of Spa2p by Pea2p variants in the El 278b background. For this analysis, 
chromosomal PEA2 and pea2-IVI409L were tagged at the carboxy terminus by integration of a cassette encoding ProA. Spa2p tagged at its amino 
terminus with the hemagglutinin (HA) epitope was expressed from a plasmid using a copper-inducible promoter. The binding affinity between Pea2p 
and Spa2p was revealed by IgG pull-down using Pea2p-ProA as bait. Pea2p and Spa2p levels were normalized against levels observed in total cell 
lysate. The ratio of recovered Spa2p by the S288C-encoded Pea2p-L409 variant is presented relative to the ratio of Spa2p recovered by El 278b- 
encoded Pea2p-M409; thus, the boxed ratios indicate that Pea2p-L409 recovered 60% the amount of Spa2p co-immunoprecipitated with Pea2p- 
IVI409. 

doi:1 0.1 371 /journal.pgen.1 004570.g003 



To assess the impact of this allelic variation on mitochondrial 
fiinction, we grew the SKI mutant with the BY4741 allele of 
MDM32 on medium containing non-fermentable glycerol as the 
sole carbon source. As shown in Fig. 5C, the allele-swapped SKI 
mutant grows poorly in glycerol-containing media, indicating that 
oxidative phosphorylation is impaired. The structure of the 
mitochondrial network is also perturbed upon introduction of 
the BY4741 aUele oi MDM32 in the SKI background. Using the 
mitochondrion-specific MitoTracker fluorescent dye, which dif- 
fuses passively across the plasma membrane and concentrates in 
active mitochondria by membrane potential, we can visualize a 
compact and collapsed mitochondrial network in SKI cells 



containing the BY4741 -encoded Mdm32p-C182/L262 variant, 
similar to that observed in mdm32S. (Fig. 5D). MDM32 is a 
paralog oiMDM3 1 , and the encoded proteins have been found to 
interact, albeit transiently and weakly, as components of protein 
complexes at the mitochondrial inner membrane [52]. We, 
therefore, assessed the effect of allelic variation at MDM32 on 
Mdm31p binding; however, we observed no difference in the 
recovery of MdmSlp by co-immunoprecipitation between the 
respective Mdm32p variants (Figure S3). In sum, MDM32 is a 
determinant of invasive growth, and replacement of the native 
SKI allele of MDM32 with the BY4741-encoded aUele yields a 
mutant filamentous growth phenotype. 



PLOS Genetics | www.plosgenetics.org 



5 



August 2014 | Volume 10 | Issue 8 | e1004570 



Genetic Variation Underlying Invasive Growth 



BY4741 SKI 

©X© 



mate & 
sporulate 



screened 1 00 spores 
from tetrads showing 
complete meiosis 



51 spores with strong invasive phenotype (++) 
24 spores with non-invasive phenotypes ( - ) 



SIA1 (581864), RUP1 (584519-586077), VORJaSC (587197-587207), SFL1 (587207-589133), /4RPS (590261-592386), ELG1 (603167-604805), 
MD/W32 (6071 52-607393), SPP2(608819), SMP3 (609994-610497), MRPL23 (6^2^24), PDR5 (619862), SLP1 (625879) 



Q 

g 4 

-' 0 



J Chr. XV 































1 100,000 200,000 300,000 400,000 500,000 600,000 700,000 800,000 900,000 1,000,000 

chromosomal position (nt) 



10 



z ) Chr. XV 



MDM32 



o 



T *. 

r 



fT+-5 

T' J* 

•.I 

I 



548,000 626,000 

Chr. XV chromosomal position (nt) 



MDM32 locus 
with allelic 
differences 
exhibiting high 
LOD scores: 



602,000 
10 



8 5 



0 



-(9.0)- 



613,000 
10 



(5.7) 




<5-6)(4.9) (S-^) 



1 



I 



Chr. XV 



MDM32, 



MR^L23 



ELG1, MDM32, SPP2, 
SMP3, MRPL23 



ELG1 



SPP2 SMP3 







BY4741 position 


nt allele change 


amino acid change 


gene 


chromosome 


(nt) 


(BY4741 to SKI) 


{BY4741 to SKI ) LOD score 



SFL1 Chr. XV 587207,588023,588410, G-A, G-A, C-T, E-L, V-M, Q-stop, 7.8 

588537,589133 T-C, A-G L-P, T-A 

MDM32 Chr. XV 607152,607393 G-C, A-T C-S, L-F 9.0 



Figure 4. Identification of allelic linkage with invasive growth from the BY4741-by-SK1 cross. A) Summary of phenotypic screening 
results. B) A chromosomal plot indicating non-synonymous allelic variation and LOD scores is presented for chromosome XV, highlighting a region of 
clustered high LOD scores (>4) around the SFLl and MDM32 loci. This region is boxed in red, and the constituent genes are indicated with the 
nucleotide positions of the allelic difference; for genes that contain multiple allelic differences, a nucleotide range is listed encompassing the 
respective residues. C) The MDM32 locus is highlighted in a plot of LOD score by chromosomal position, and the inset box indicates the orientation of 
MDM32 and neighboring genes on chromosome XV. LOD scores are indicated for each gene in the boxed region, and as in Figure 2, the width of the 
LOD score bar indicates the chromosomal boundaries of the allelic changes. D) Specific nucleotide and amino acid changes between BY4741- and 
SKI-encoded alleles of the SFLl and MDM32 genes are indicated, along with the respective LOD score for each gene. 
doi:1 0.1 371 /journal.pgen.1 004570.g004 



Discussion 

The linkage analysis presented here identifies a broad gene set 
contributing to pseudohyphal growth (Fig. 6). The pattern of 
allelic linkage indicates a large number of determinant loci 
underlying invasive growth, consistent with results from systematic 
single-gene deletion and overexpression studies [28,31,32]. Within 
this gene set, components of the polarisome and mitochondria play 
important roles in enabling invasive growth. We report here that 
Pea2p residue 409 impacts bud site selection in haploid cells 
undergoing invasive growth, although the effect is less pronounced 
in determining initial distal-versus-proximal budding in virgin 
mother cells. Pea2p residue 409 lies distinct from the Pea2p coUed- 
coil region between residues 236 and 327 and is important for 
Spa2p binding. Spa2p interacts with Stellp and Ste7p from the 
Ksslp MAPK pathway, providing a mechanism for polarisome- 



mediated regulation of Ksslp MAPK activity [41]. Our results 
further indicate that nuclear-encoded Mdm32p is required for 
invasive growth in SKI, and that residues 182 and 262, located 
outside of the mitochondrial pre-sequence (AA 1-102) and at the 
boundary or outside of a transmembrane domain (AA 161-184 
and 636-653), are important in enabling invasion, as well as in 
determining aerobic respiratory function and mitochondrial 
morphology. Mdm32p is proposed to function cooperatively with 
other inner membrane proteins and components of the outer 
mitochondrial membrane in the maintenance of mitochondrial 
morphology [52], potentially through cytoskeletal interactions that 
may be affected by variation at these sites. 

This analysis highlights two additional points. First, the non- 
filamentous BY4741 background is not uniformly repressive with 
respect to pseudohyphal growth. From the SI 278b and SKI 
crosses, we identified three and five blocks of allelic variation. 
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Figure 5. Phenotypic consequences of variation in Mdm32p residues 182 and 262. A) Invasive growth assays of wild-type and mdm32 
mutants in the SKI background were performed as described previously. Invasive growth levels were quantified by determining the ratio of the pixel 
intensity of images of the spotted cultures before and after plate washing, and these ratios are indicated in the boxes to the lower right of the 
spotted cultures post-washing. B) Introduction of the BY4741 -encoded MDM32-C1 82/L262 allele into SKI results in decreased cell elongation relative 
to wild type. In this analysis, cells were scraped from colonies into a suspension for imaging of cell morphology as described previously. C) A strain of 
SKI carrying the S288C-encoded allele of MDM32 exhibits a growth defect on medium containing the non-fermentable carbon source glycerol. Cells 
from the indicated strains were grown overnight in glucose- containing medium, and subsequently, 10-fold serial dilutions were spotted onto plates 
containing glucose (YPD) or glycerol (YPG) as a carbon source. YPD plates were incubated for 2 days and YPG plates for 3 days at the indicated 
temperatures. D) Analysis of mitochondrial network structure in a mutant strain of the SKI background carrying the BY4741 -encoded MDM32-C182/ 
L262 allele. As visualized by staining with a mitochondrial dye, the mitochondrial network is compact and collapsed in a strain deleted for MDM32 and 
in the SKI mutant carrying the BY4741 allele. 
doi:1 0.1 371 /journal.pgen.1 004570.g005 



respectively, in BY4741 linked with the invasive growth pheno- 
type; full listings of the encompassed alleles with respect to each 
cross are presented in Tables S5 and S6. Alleles in S288C-derived 
strains that promote pseudohyphal growth are antagonized by 
alleles that repress filament formation, such as the pseudogene 
form of FLOS; similarly, alleles in SKI linked with the non- 
invasive phenotype may be offset by alleles that promote invasion, 
such as the SFLl allele containing a premature stop codon. 
Second, the identified allelic variation in PEA2 and MDM32 and 
the aUele-swapping experiments performed here indicate that 
within a given genome, functionally interacting genes coevolve to 
impact phenotype. The majority of genetic variation linked with 
invasive phenotype in this study involves site-specific changes that 
do not create pseudogenes. Alleles of these genes yield functional 
proteins within the respective genomic contexts; however, a given 
allele results in a hypomorphic phenotype when introduced into a 
non-native strain. It should be noted that the BY4741 -encoded 
allele of PEA2 may be viewed as being pseudohyphal competent, 
as Liu et al. [36] reported that the introduction of Z 1278b- 
encoded FLOS in a S288C-derived strain is sufficient to enable at 
least some degree of pseudohyphal growth. The data here suggest 
tiiat partner genes have likely co-evolved with genes such as PEA2 
and MDM32, and the resulting protein complexes are, thus, an 



important determinant of cell phenotype. These findings highlight 
the utiUty in studying these complexes as a whole, in supplement to 
individual proteins, in order to accurately identify the functions 
and properties that specify phenotype. 

It is interesting that the studies here indicated very little overlap 
between alleles linked with invasive growth in the S1278b and 
SKI strains with respect to BY4741 (Fig. 6). Studies mapping 
quantitative trait loci (QTL) in a cross of the laboratory strain 
BY4716 and the vineyard strain RM-11 identified hotspots 
impacting gene expression, protein abundance, and small mole- 
cule-dependence [55-58] . These hotspots were principally due to 
alleles in the BY4716 background, leading Ronald and Akey [59] 
to suggest that the causative polymorphisms may occur at low 
frequency. The non-overlapping allele sets identified in our 
analysis are not suggestive of hotspots, but rather highlights the 
substantial importance of epistatic interactions in determining the 
sum filamentous phenotype resulting from variant alleles in the 
haploid segregants. These epistatic interactions likely represent 
instances of gene coevolution, which has been suggested to occur 
at an elevated rate for genes encoding proteins of shared biological 
functions and/ or for proteins that have coevolved between species 
[60,61]. Clark et al. [60] have analyzed the rate of covariation for 
pairs of proteins over evolutionary time, and by this analysis. 
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doi:1 0.1 371 /journal.pgen.1 004570.g006 



polarisome components as a whole do not exhibit statistically 
significant evidence of covariation, although many mitochondrial 
complexes do yield signature indicating evolutionary rate covari- 
ation. Further analyses of individual protein pairs from the strains 
used in our study will be necessary to identify a set of coevolved 
proteins that drive the filamentous growth phenotype. 

In summary, we used pooled segregant whole-genome sequenc- 
ing to dissect gene networks that determine yeast pseudohyphal 
growth. This analysis identified allelic variation in the known 
pseudohyphal growth genes FLOS and FLOll, while also 
revealing variation in the negative regulator SFLl, the coding 
sequence of which contains a premature stop codon in the invasive 
SKI background. We further found that amino acid 409 in the 
polarisome protein Pea2p is a site of allelic variation critical for the 
protein's ability to signal through the Ksslp MAPK pathway, 
estabhsh unipolar budding during pseudohyphal growth, and bind 



the Spa2p polarisome scaffold. Linkage analysis identifies variation 
in MDM32 as a determinant of invasive growth between S288C 
derivatives and the SKI strain; the 182 and 262 residues are sites 
of variation and contribute to Mdm32 function in aerobic 
respiration and invasive growth. 

Materials and Methods 

Yeast Strains and Growth Conditions 

A listing of yeast strains and plasmids used in this study is 
provided in Tables S7 and S8. Haploid deletion mutants were 
constructed by PCR-mediated gene disruption using pFA6a- 
KanMXG or pUG72 [62,63]. Yeast strains were propagated on 
rich YPD medium (1% yeast extract, 2% polypeptone and 2% 
glucose) medium or synthetic medium as described [64]. Yeast 
invasive growth was assayed on YPD medium. 
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Pooled Segregant Analysis 

The statistical modeling used to derive the probabilities of 
idc-ntifying linkage are described by Birkeland e.l al. [65]. 
Following mating as indicated in Fig. 1, resulting strains were 
sporulated and asci were dissected. The dissected spores were 
grown overnight at SCG and were individually tested for mating 
type. Spores resulting from complete meiosis (four viable spores 
with two each of the a and a mating types) were then used for 
whole genome sequencing. Each spore was assigned to invasive 
or non-invasive pools based on its invasive growth phenotype. 
To ensure equal representation of all segregants in a pooled 
population, each haploid strain was grown overnight at 30°C in 
individual 4 ml YPD cultures. The ODeoo of the cultures was 
determined and used to calculate the appropriate volume of each 
strain so that upon mixing, we would achieve equal numbers of 
cells. A mate-pair library with 300-bp fragments was prepared 
for each of the phenotypic pools, and each library sequenced as 
paired-end reads using the lUumina Genome Analyzer (Univer- 
sity of Michigan DNA Sequencing Core). Sequence analysis was 
performed as described [65]. To obtain an estimate of the 
recombinant and non-recombinant spore counts in each 
phenotypic pool for a given observed sequence variant, we 
multiplied the number of spores in the pool by the fraction of 
sequence reads from that pool that matched the corresponding 
allele variant. These values were then used in standard LOD 
score calculations. 

Analysis of Yeast Budding Pattern 

Budding patterns of invasive cells were determined as 
previously described [44,45]. In brief, equal concentrations of 
mid-log phase cultures were spotted onto YPD plates and 
incubated for 7 days at 30°C; surface cells were subsequently 
washed off under a gentle stream of water. Residual invaded cells 
were recovered from the agar using a sterile toothpick. Cells were 
washed twice in sterile water and were stained with 2 ng/ml 
calcoflour white. Bud scars were visualized by fluorescence 
microscopy. Cells with more than three bud scars were examined. 
Budding patterns were determined by criteria previously 
described [45]. Budding patterns were divided into four sub- 
groups: axial, bipolar, unipolar and random. The axial pattern 
was defined as a long chain of bud scars on the proximal cell pole. 
Cells with a cluster of scars exclusively at the distal pole were 
classified as exhibiting unipolar budding. A pattern of medial bud 
scars was scored as random budding, whereas cells with bud scars 
erjually distributed on both proximal and distal poles were 
classified as undergoing bipolar budding. For these analyses, 200- 
250 cells from each strain were scored. 

Plasmid constructioii. For expression of S1278b alleles in 
nuU mutants, yeast open reading frames of the Y826 background 
with 1 kb of upstream sequence and 300 bp of downstream 
sequence were cloned into a Gateway vector by standard methods 
of recombination-based Gateway cloning (Invitrogen) [66,67]. For 
fusing enhanced green fluorescent protein (eGFP) and tandem 
affinity purification (TAP) cassettes to the carboxyl-terminus of the 
budding yeast Mdm3 Ip and Mdm32p proteins, Gateway plasmids 
416-GPD-ccdB-EGFP and 414-GPD- ccdB-TAP were used [68]. 
Plasmid pCu-Spa2p-3xHA-[/i?/15 was modified from pCu-3xHA- 
URA3 using standard restriction enzyme digestion and ligation- 
based techniques. The Sail and Hindlll restriction sites of this 
plasmid was used to integrate the SPA2 open reading frame 
between the copper-inducible CUPl promoter and 3xHA/GFP 
tag. 

Integrated point mutations. Allelic variants of PEA2 were 
generated in the SI 278b background as integrated point mutations. 



A URA3 cassette was amplified from the pUG6 plasmid and then 
used to replace the PEA2 open reading frame in Y826. 
Subsequentiy, 5-FOA-mediated counter selection [69] was applied 
to replace the integrated URA3 cassette with a 500 bp cassette 
encompassing the PEA2 allelic change (M409L), recreating the 
S288C-derived allele in the SI 278b background. The integrated 
point mutation was confirmed by sequencing of an amplified PCR 
fragment using the University of Michigan Sequencing Core. 

Yeast invasive growth assays. Invasive growth of haploid 
strains was determined by the standard plate washing assay of 
Gimeno et al. [1]. Mid-log phase cultures were spotted onto YPD 
plates and incubated for 7 days at 30°C, and surface cells were 
subsequently washed off under a gentle stream of water. Residual 
cells from the spotted culture on the plate were imaged using a 
Nikon Eclipse 80i upright fluorescence microscope. 

FKE-lacZ assays. Plasmids pLG669-Z FLOl 1-6/7 and -9/ 
10 were transformed into designated yeast strains by standard 
protocols [64] . Cells were grown in SC-Ura media overnight, and 
then inoculated into fresh media to an ODgoo of approximately 
0.2. Cultures were subsequently grown for 3-4 hours to an ODgoo 
of approximately 0.8-1.0. P-galactosidase activity was determined 
with the use of ortho-nitrophenyl-fi-galactoside as a substrate 
(Sigma Aldrich) [70,71]. 

Assays for respiratory activity. Respiratory ability was 
assessed using YPG plates with glycerol as a non- fermentable 
carbon source. Cells were pre-grown overnight and then inoculated 
into fresh SC medium to an ODhoq of approximately 0.2; culture 
were subsequently incubated an additional 3 hours to an ODroo of 
approximately 0.8. The same amount of cells were serially diluted 
and spotted onto YPD and YPG (1 % yeast extract, 2% polypeptone 
and 2% glycerol) plates, and incubated for 2 days at 30°C. 

Fluorescence microscopy and image processing. For 
live-cell imaging, spotted cultures were grown in YPD plates for 
7 days at 30°C. After washing the surface of the plates with a 
gentle stream of water, cells that had invaded the agar were 
recovered using a steel scalpel, and the extracted cells were washed 
twice with distilled water. Images were taken using a Nikon Eclipse 
80i upright fluorescence microscope. To determine budding 
patterns, cells were stained with Calcofluor White as previously 
described [45]. Over 250 cells were counted twice for each strain. 
Statistics were analyzed using the ANOVA approach. To reveal 
mitochondrial morphology, living cells were stained vnth 50 nM 
MitoTracker Red, and the staining pattern was visualized by 
fluorescence microscopy [72]. 

Immunoprecipitation studies. PEA2 alleles were generat- 
ed by tagging with ProA on the C-terminus in wild type and pea2- 
M409L mutants in the SI 278b strain. The SPA2 open reading 
frame was cloned into pCu-3HA-URA, such that SPA2 was 
expressed from a copper-inducible promoter as a fusion protein 
with three copies of the hemagglutinin epitope at its amino 
terminus. The resulting plasmids pCu-3HA-ScSPA2 and pCu- 
3HA-SSPA2 were introduced by transformation into both wild 
type and pea2-M409L mutants. The binding aflinity between 
selected proteins was detected by co-immunoprecipitation. For 
native immunoprecipitation, 2 OD units of cells were lysed in 1 ml 
lysis buffer (50 mM Tris-HCl [pH 7.5], 150 mM NaCl, 2 mM 
EDTA, 0.5% Triton X-100, 1 mM PMSF, and Complete EDTA- 
free protease inhibitor [Roche]) with glass beads. After centrifu- 
gation at 13,000 g for 10 min, the resulting supernatant was 
incubated with protein G-Sepharose 4 Fast Flow (GE Life Tech) 
for 2 hr at 4°C. After washing the Sepharose with lysis buffer six 
times, the bound materials were eluted by boiling the Sepharose in 
SDS-PAGE loading buffer. The resulting eluate was analyzed by 
Western blotting with designated antibodies (anti-HA, anti-Protein 
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A and anti-GFP). Blots were developed using the SuperSignal 
West Dura Extended Duration Substrate (Thermo Scientific). 

Supporting Information 

Figure SI Functional characterization of genes exhibiting allelic 
linkage with invasive growth in the BY4741 cross with 21278b. 
Genes from the analysis that have been previously identified as 
exhibiting filamentous growth phenotypes upon deletion [32] are 
listed to the left. The remaining genes are listed to the right, with 
genes of unknown functions and genes with functions that may be 
related to pseudoh^phal growth indicated separately. Functions 
are drawn from data deposited in the Saccharomyces Genome 
Database as of publication. The chromosome in which the gene is 
located is indicated in parentheses. It should be noted that this list 
encompasses genes that may only exhibit linkage because of their 
close proximity to an important allelic determinant. 
(TIF) 

Figure S2 Functional characterization of genes exhibiting allelic 
linkage with invasive growth in the BY4741 cross with SKI. Genes 
from the analysis that were identified as exhibiting filamentous 
growth phenotypes upon deletion [32] are listed to the left. The 
remaining genes are listed to the right, with genes of unknown 
functions and genes with functions that may be related to 
pseudohyphal growth indicated separately. Functions are drawn 
from data deposited in the Saccharomyces Genome Database as of 
publication. The chromosome in which the gene is located is 
indicated in parentheses. It should be noted that this list 
encompasses genes that may only exhibit linkage because of their 
close proximity to an important allelic determinant. 
(TIF) 

Figure S3 Binding of MdmSlp by BY4741- and SKl-encoded 
Mdm32p variants. The MDM32 open reading frames from SKI 
and BY4741 were cloned into a GPD-eGFP Gateway plasmid, 
such that the MDM32 sequence was expressed from the GPDl 
promoter as an in-frame 3'-fiision to sequence encoding enhanced 
GFP. MDM31 was cloned into the GPD-TAP plasmid, yielding a 
fusion of the tandem affinity purification (TAP) tag to the carboxy 
terminus of MdmSlj) upon (expression from the GPDl promoter. 
The resulting plasmids were transformed into mdm32A mutants. 
The binding affinity between Mdm32p-GFP and Mdm31p-TAP 
was revealed by IgG pull-down using Mdm31p-TAP as bait. No 
significant difference in binding between Mdm31p-TAP and the 
respective Mdm32p-GFP variants was observed. 
(TIF) 
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Table SI Genes exhibiting allelic linkage (LOD>3) from the 
BY4741-by-21278b cross. AUeles in a single chromosome are 
demarcated by a double line; groupings of alleles within a possible 
linkage block are separated by a single line. 
(DOCX) 
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